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Abstract 

The Earth effects on the energy spectra of supernova neutrinos are studied. We 
analyze numerically the time-integrated energy spectra of neutrino in a mantle- 
core-mantle step function model of the Earth's matter density profile. We consider 
a realistic frame- work in which there are three active neutrinos whose mass squared 
differences and mixings are constrained by the present understanding of solar and 
atmospheric neutrinos. We find that the energy spectra change for some allowed 
mixing parameters. Especially, the expected number of events at SNO shows char- 
acteristic behavior with respect to energy, i.e. a great dip and peak. We show that 
observations of the Earth effect allow us to identify the solar neutrino solution and 
to probe the mixing angle 92- 
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1 Introduction 



Neutrino oscillation is widely accepted as a solution to the solar and atmospheric neutrino 
problems. However, in spite of the remarkable progress made in the studies of solar and 
atmospheric neutrinos, there are some ambiguities about neutrino mass spectrum and 
mixing. Then neutrino emission from collapse-driven supernova attracts a significant at- 
tention because the supernova neutrinos can probe whole neutrino mass and their studies 
will help to identify the pattern of neutrino mass and mixing. In fact, the direct ob- 
servations of neutrinos from SN1987A by the Kamiokande and the 1MB detectors 
have provided rich information not only on the mechanism of supernova but also on the 
neutrino physics 0, Q and paved the way for a new phase in the neutrino astrophysics. 

i 

The neutrinos, which are produced in the high dense region of the iron core, ^ interact 
with matter before emerging from the supernova. The presence of non-zero masses and 
mixing in vacuum among various neutrino flavors results in strong matter dependent 
effects, including conversion from one flavor to another. Hence, the observed neutrino 
flux in the detectors may be dramatically different for certain neutrino flavors and for 
certain values of mixing parameters, due to neutrino oscillation. ^ ^ 

Similar matter effects occur in the Earth, too (in connection with solar neutrino and 



SN1987A, see |T0[ and [|TT|, |T2|, |T3|, respectively). The neutrino trajectory inside 
the Earth before reaching the detector depends on the direction of the supernova relative 
to the Earth and the time of the day. The comparison of signals from different detectors 
would allow one to reveal the Earth matter effects. Also certain features of the energy 
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spectra can reveal the Earth matter effect even from the observations in one detector. In 
this paper, we study numerically the Earth matter effects on supernova neutrino spectra 
with all three active neutrinos considered. We show that some ambiguities in the solar 
neutrino solutions and in 62 can be resolved by the Earth effect of supernova neutrino. 

2 Three-flavor formulation 

We consider three-flavor mixing among active neutrinos. The three flavor eigenstates are 
related to the three mass eigenstates in vacuum through a unitary transformation, 

3 

\^a) = Uaa\l^a), (1) 
a=l 

where |i/q,) are flavor eigenstates with flavor a{a = e,fi,T) and |z/a) are mass eigenstates 
with mass m^. A convenient parameterization for U = ^2) ^3) is given by 



U 



^ C2C3 S3C2 S2 

~S3Ci — S1S2C3 C1C3 — S1S2S3 S1C2 I , (2) 

V S1S3 - S2C1C3 -S1C3 - S2S3C1 C1C2 



where Si = sinOi and q = cos6'j for i = 1,2, 3. The quantities 6i are the vacuum mixing 
angles. We have here put the CP phase equal to zero in the CKM matrix. 

Experimental evidence suggests [1^, 0, [l^ that both solar and atmospheric neutrinos 
display flavor oscillations, and hence that neutrinos have mass. All the existing experi- 
mental results on the atmospheric neutrinos can be well described in terms of the z/^ Uj- 
vacuum oscillation with mass squared difference and the mixing angle given by 

^^ata ~ 7 X lO^^eV^ , sin^ 29 ^ 1. (3) 
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In contrast, the current (mean event rate) solar neutrino data admit three types of 
MSW solutions if the solar v^. undergo two-neutrino transitions into active neutrinos, 



^ [0: the SMA, the LMA and the so-called "LOW" solution. We use SMA and 



LMA parameters as samples. 

(SMA) Am^ ^ (4 ~ 10) x lO^'^el/^ , ^m^^e^ ^ (2 ~ 10) x 10"^ (4) 
(LMA) Am| ^ (1 ~ 10) x 10~^e1/2 , sin2 2^0 ^ 0.7 ~ 0.95 (5) 

These ambiguities mentioned above can be resolved by supernova data including the Earth 
effects, ra 



3 Calculating the Earth effects 

The calculations of the energy spectra of neutrinos performed in this paper are based 
on the time-evolution-operator method developed by Ohlsson and Snellman. This 
method provides the time evolution operator in terms of the mass squared differences and 
the vacuum mixing angle without introducing the auxiliary matter mixing angles. The 
evolution operator Uf{L, A, E) is a 3 x 3 matrix whose components are functions of the 
length L of the neutrino trajectory, the matter density parameter A and the neutrino 
energy E. The neutrino conversion probability amplitude is given by 

A,^^,^ = {iyf,\UfiL,A,E)\u^), (6) 

where A = Gpp/ V^rnN, a, P = e, /i, r, Gp is the Fermi constant, m^v is the nucleon mass, 
and p is the matter density. 
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When the neutrinos travel through a series of matter densities with matter density 
parameters Ai, A2, ■ ■ ■ , An and thicknesses Li, L2, ■ ■ ■ , Ln, the total evolution operator is 
simply given by 

Uf" = Uf{Ln, An, E) ■■■ Uf{L2, A2, E) Uf{L,, A,, E) (7) 

We use the following values as the matter density parameter of the Earth core and mantle, 
respectively; 

A,ore = 4.35 X , {p^re = ILSg/cm^), (8) 

Amantle = 1-70 X IQ-^^cV , {pmantle = 4.5g/cm2). (9) 

When neutrinos from such a distance as in the case of supernova arrive at the Earth, 
they are not coherent and interference between the different mass eigenstates does not 
occur 1^ . This effect can be included in our calculation by time-averaging the conversion 
probabilities at the surface of the Earth. 

In order to compute the neutrino energy spectra, we used energy spectra of neutrinos 
arriving at the surface of the Earth (Fig.|l|) calculated by Watanabe et al. [^] They have 
been studying the future detection of a supernova neutrino burst taking the three-flavor 
neutrino oscillation inside the star into account. The system has two resonances: one at 
higher density (H- resonance) and the other at lower resonance (L-resonance). Adiabaticity 
at these resonances significantly influences on not only the neutrino flux outside the star 
but also the magnitude of the Earth effects. 
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We computed energy spectra of Ve and for the foUowing sets of mixing parameters, 
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where mi < < (normal hierarchy) and Am^j = "^2 ~ ""^L ^^23 ~ "^3 ~ ^2- '^^^ 
value of 62 determines the adiabaticity at the H resonance. [l[] The conversion is adiabatic 
and nonadiabatic for sin^ 202 = 0.04 and 1 x 10"^ respectively, sin^ 2^3 = 0.8 and 7 x 10-3 
correspond to the solutions to the solar neutrino problem LMA and SMA, respectively. 

In terms of the classification by Dighe and Smirnov [|T|, the parameters (a), (b), (c) 
and (d) correspond to I-LMA, I-SMA, III-LMA and III-SMA, respectively. They predict 
that the Earth effect would be observed in spectra of z/g of (a), both u^. and of (c) and 
i>e of (d). In our calculation, on the contrary, significant Earth effects (> 10%) can be 
seen in three among 2x4 spectra stated above: of (a) and both and z/g of (c), not 
i/g of (d). 

The flux reaching the detector from a supernova at a distance d from the Earth is 
reduced by an overall geometric factor of l/(47r(i^). We used 10 kpc as the value of d. 
We computed the flux of neutrinos at 130 discrete energies in a mantle-core-mantle step 
function model of the Earth's matter density profile. For each neutrino flux with specific 
energy, we performed a sequence of energy averaging operations as follows. First we 
compute the z/q —>■ conversion probabilities Py^^yi^{E) at 0.1 MeV intervals. Then the 
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flux at E MeV at the detector F^{E) is, 

E+lMeV 

i^.^(^)= E [^^. E ^..^..(^')], (10) 

flavcyr E'=E-lMeY 

where Fj,^ is the z/q flux just arriving at the Earth. We take the energy width 2 MeV 
taking a rough energy resolutions of detectors into account. In fact, the energy resolutions 
of detectors depend on energy it self. But the rough estimation of energy resolution is 
sufficient, since our purpose is to see overall shapes of the spectra. 

4 Results and discussions 

We calculated energy spectra of neutrino which traveled through the mantle(3000 km)- 
core(7000 km) -mantle(3000 km), which correspond to zero nadir angle. Significant Earth 
effects were seen in the flux of with mixing parameter (c) and z/g with mixing parameter 
(a) and (c). Then we converted the fluxes in which the Earth effects can be seen to energy 
spectra at SNO and SuperKamiokande. In calculating the energy spectra at SNO, two 
charged current interactions, 

Ve + d p + p + e~ (11) 
Ue + d n + n + e^ (12) 



are taken into account. For the cross sections for these reactions, we refer to p2|. The 
detection efficiency is taken to be one. Fig.0 shows the time-integrated energy spectra of 
I'eyi^e sud z/g + at SNO with mixing parameters (c) with and without the Earth effects. 
In calculating the energy spectra at SuperKamiokande, the reaction, 

iye+ P ^ n + e~^, (13) 



which dominates the total event number, is taken into account. For the cross section of 
this reaction, we refer to The appropriate detection efficiency curve is also taken 

into account |2^. Fig.^ shows the time-integrated energy spectra at SuperKamiokande 
with mixing parameters (a) and (c) with and without the Earth effects. We don't show 
the energy spectra at SNO with mixing parameter (a) because the Earth effects on z/g 
and Ve are negligible and small, respectively, and the number of event coming from is 
dominant. For mixing parameter (b) and (d), no significant Earth effect can be seen. 

The flux of the neutrino passing through the Earth oscillates with respect to the energy. 
This behavior can be explained as follows. For simplicity, we assume the density of the 
Earth be homogeneous. The dashed line in Fig.^ shows the flux of Ve which traveled 4000 
km through the Earth core. Probability of neutrino conversion, for example, 
can be expressed as a function of the length the neutrino trajectory in matter and the 
neutrino energy, 

P.,_^ = sin^2^Msin2y^, (14) 

where Lm{E) is the matter oscillation length which can be written with the vacuum 
oscillation length Ly, vacuum mixing angle 6v, and the neutrino-electron interaction 
length Le as 

^''^^^ = {i-2^cos2ey + ii^yy/^' ^^^^ 

where Ly = A'kETi/ ^m?c^ , Le = \/2nhc/GFne, and Ue is the electron number density. 
Substituting, 

R = 4000km , cos20v' = , rig = Ucore, (16) 
we obtain the argument of the second factor of Py^^y^, i6Vi+ao22 5£g.^^ This function 
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decreases from 13 to 8 as increases from 5 MeV to 15 MeV and takes an almost 
constant value 8 for E > 15 MeV: the oscillation length asymptotically approaches a 
constant value. (Fig. |) 

On the other hand, the difference between the flux at the detector and the u^, 
flux Fe before traveling the Earth is proportional to, [Q 

- oc FO - F°, (17) 

where Fg" and F° are the original fluxes of, respectively, and z/^,,- at the center of the 
supernova, which vanish at F 20 MeV. |21|] Fig.|l] shows the original fluxes: the solid, 
dashed and dotted lines show Fg°, Fgand F°, respectively. 

The shapes of spectra can be understood roughly by the latter factor. At lower energy 
FO - F" > 0, then Ff - F^ > 0. F° - F° = at F ^ 20 eV , then Ff - F^ = 0. At 
higher energy Fg° — F° < 0, then F^ — Fg < 0. The small oscillation can be understood 
as the oscillation of P^^^^^. 

Fig.|^ shows the z/g flux inside the Earth with parameters (c) which travel the mantle, 
core, and mantle. The length of the mantle and core are set to be 3000 km and 7000 km, 
respectively. If averaging operations are performed, flux oscillations damp because u flux 
of different energies have different oscillation lengths and cancel the oscillation behavior 
while traveling long distance. Since the oscillation lengths are comparable to the radius 
of the Earth, the Earth effects are strongly dependent on the neutrino path inside the 
Earth. So the magnitude of the Earth effects will be different at each detector. 

Significant Earth effects can be seen in the energy spectra at SNO for mixing parameter 
(c) (Fig.^. There are a large dip at 30 MeV and a peak at 37 MeV. Binning appropriately, 
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(for example 25 ~ 35 MeV and 35 ~ 40 MeV, etc.) we can tell whether the Earth effects 
exist. Since the energy value of dip and peak change if zenith angle of neutrino path 
changes, appropriate binning will change accordingly. 

The Earth effects in Fig.§ are very small. This smallness comes partly from the small- 
ness of Fg — (Fig.|l]). Although the expected number of events at SuperKamiokande 
is very large (~ 10000 counts) |2l| and then statistical error would be small, it would be 
difficult to detect the Earth effect on neutrinos from supernova at lOkpc from the Earth 
due to the uncertainties in supernova model. In this case we suggest as a criterion of 
the Earth effect the ratio of events at lower energy (< 20 MeV) to that at higher energy 
(> 20 MeV). These subjects need quantitative and statistical argument, which we plan 



to study in great detail elsewhere. p6 



If we observe a supernova and detect neutrinos in future, whether the Earth effects 
exist can help to determine the neutrino parameters. 

• When the Earth effects are observed at both SuperKamiokande and SNO, the neu- 
trino parameter is (c). 

• When the Earth effects are observed at SuperKamiokande only, the neutrino pa- 
rameter is (a). 

Recent result of solar neutrino analysis at SuperKamiokande favor LMA rather than 



SMA. If the solution of the solar neutrino problem would turn out to be LMA, the 
Earth effects on supernova neutrino can distinguish parameter (a) from (c): we could 
determine 62, which is only loosely constrained by the other current experiment. 
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Figure 1: Energy spectra of i^e, Ve and of the numerical supernova model used in this 
paper. The solid, dashed and dotted lines are the spectrum of i^e, and v^i respectively. 
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Figure 2: Time-integrated energy spectra at SNO with parameter (c). Short-dashed and 
dotted hne correspond to events without and with the Earth effects, respectively. Long- 
dashed-dotted and short-dashed-dotted hne correspond to Pg events without and with the 
Earth effects, respectively. Solid and long-dashed line correspond to total events without 
and with the Earth effects, respectively. 
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Figure 3: Time-integrated energy spectra at SK. Solid line corresponds to energy spectra 
without the Earth effect. Dashed and dotted line correspond to energy spectra with Earth 
effect using parameter (a) and (c), respectively. 
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Figure 4: Energy spectrum of v^. with parameters (c) which did not travel (sohd) and 
traveled through the core 4000km (dashed). 
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Figure 5: Argument of the second factor of the probabihty P^^ 
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Figure 6: flux in the Earth. The sohd and dashed hnes show z/g flux at 10 MeV 
without and with the averaging operations, respectively. The dotted and dash-dot-dash 
hnes show flux at 30 MeV without and with the averaging operations, respectively 
(almost degenerate). 
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